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Thin-film sub-5 nm magnetic skyrmions constitute an ultimate scaling alternative for future digi¬ 
tal data storage. Skyrmions are robust non-collinear spin-textures that can be moved and manip¬ 
ulated by small electrical currents. We show here an innovative technique to detect isolated nano- 
skyrmions with a current-perpendicular-to-plane geometry, which has immediate implications for 
device concepts. We explore the physics behind such a mechanism by studying the atomistic 
electronic structure of the magnetic quasiparticles. We investigate how the isolated skyrmion 
local-density-of-states which tunnels into the vacuum, when compared to the ferromagnetic back¬ 
ground, is modified by the site-dependent spin-mixing of electronic states with different relative 
canting angles. Local transport properties are sensitive to this effect, as we report an atomistic 
conductance anisotropy of up to ~20% for magnetic skyrmions in Pd/Fe/lr(111) thin-films. In sin¬ 
gle skyrmions, engineering this spin-mixing magnetoresistance possibly could be incorporated in 
future magnetic storage technologies. 


1 



Si complimentary-metal-oxide-semiconductor 1 (CMOS) compatible magnetic devices rep¬ 
resent the current state-of-the-art in information data storage circuits 2 . In such devices, the 
information is encoded by manipulation of different spacial magnetic domains, and the data 
is read by sensing the variation in the magnetoresistance as a function of the magnetization 
direction 3 . 

Incorporation of magnetic skyrmions instead of domain walls (DWs) may improve de¬ 
vice performance and scaling possibilities 4 . DWs are sensitive to defect-pinning 5- ' while 
skyrmions intrinsically are not 8-10 . This insensitivity to defects explains why skyrmions 
can be moved by electric current densities orders of magnitude smaller than DWs 11,12 while 
achieving smaller yet still sizeable magnitudes of velocity 10,13-15 . Skyrmions, for topological 
reasons 16 , are relatively robust magnetic particle-like configurations 17,18 , stable even up to, 
so far, near-room temperature 19 . They can be confined at will and their shape and size 
controlled with an external magnetic held 20 . 

An illustrative device scaling route is an extension 4 of the racetrack memory 21 to in¬ 
corporate single 22,23 confined magnetic skyrmions instead of DWs. Such a concept con¬ 
stitutes an ultimate scaling alternative in terms of packing density, speed, and power 
consumption 4 . We envision a thin-hlm magnetic heterostructure 24 where sub-5 nm chiral 
skyrmionic quasiparticles are generated (via materials engineering and external magnetic 
fields) and moved laterally along a magnetic racetrack by in-plane currents, which has been 
shown experimentally 11,12,25 . The bit-wise data would be encoded by out-of-plane currents 
which create or annihilate the individual skyrmions, thus setting or resetting the bit-state. 
This has also been shown experimentally 23 ; thus, two important ingredients to a viable 
skyrmion racetrack memory system, lateral bit-wise movement and set /reset of each bit- 
state, have been established experimentally. 

But how could one read each bit-state? Current-in-plane (CIP) detection of skyrmions has 
been shown experimentally 25-27 and understood theoretically as a topological Hall effect 28,29 , 
but may be costly in terms of power consumption and difficult to fabricate in terms of device 
geometries. A better option would be the direct detection of the nano-skyrmions via current- 
perpendicular-to-plane (CPP) 30,31 geometry. 

We demonstrate that such a detection mechanism for single skyrmions is indeed possible, 
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Figure 1 | Perpendicular detection of nano-skrymions. (a) Illustrative concept device for 
perpendicular reading of single nano-skyrmions. Due to energy-dependent spin-mixing 
perturbations in the atomistic electronic structure as a function of position within skyrmions, 
the electric current relation i FM ^ i Sk holds, and therefore magnetic data information can be 
sensed in a CPP-geometry. (b) Illustrative STM-spectroscopy experiment of fcc-Pd/Fe 
overlayer on single-crystal fcc-lr(lll) bulk substrate. The tunneling conductance is modified by 
the combined effects of local magnetic non-collinearity and substrate-induced spin-orbit 
interaction. For similar physical reasons as (a), the tunneling conductance at position xq is 
different from that at position x\. 

even using non-spin-polarized electric currents. Using Fig. la as a schematic illustration, we 
will show how the tunneling current between a suspended metal contact through vacuum 
depends on the non-collinear magnetic state-of-phase below. (Suspended metal contacts are 
possible with state-of-the-art fabrication techniques 32 , but one could also imagine tunnel¬ 
ing through a weakly-interacting two-dimensional (2D) insulator, such as hBN or M 0 S 2 ). 
Such a process can be intimately understood in a non-spin-polarized scanning tunneling mi¬ 
croscopy (STM) experiment, Fig. lb. This all-electrical detection is a departure from typical 
experiments which rely on spin-polarized injection to detect magnetic structures, and thus 
inherently an improvement from a device application perspective. 

We will examine how the tunneling current between a non-magnetic STM tip and surface 
varies as a function of the skyrmion’s local atomic magnetization direction. Considering the 
skyrmion’s central spin-flipped atom, its local magnetic environment is inequivalent from 
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those atoms in the ferromagnetic background. Thus, one could expect already a conductance 
anisotropy in the center of the skyrmion due to the change in the electronic structure relative 
to the ferromagnetic state. One purpose of our work here is to investigate this effect. In 
addition, the heavy-metal substrate induces a large spin-orbit interaction (SOI), coupling the 
local magnetization to the real-space direction, further modifying the electronic structure 
as a function of magnetization rotation relative to the substrate plane. This tunneling 
anisotropic magnetoresistance (TAMR) 33 effect has been studied in detail but typically for 
homogeneous magnetic domains (either different ferromagnetic domains 33,34 or homogeneous 
spin-spirals 35 ). We seek to extend this theoretical basis to inhomogeneous spin-textures such 
as skyrmions. We focus specifically on single skyrmions and do not investigate networks or 
lattices of skyrmions. 

In this Article, we explain the combined (non-collinear and SOI-induced) spin-mixing 
magnetoresistance in terms of fully self-consistent calculations from first-principles, where 
we have direct access to the electronic structure of not just the magnetic heterostructure, 
but even of the states decaying into the vacuum, which are essential to the tunneling con¬ 
ductance. We find a rather large atomistic conductance anisotropy of up to ~20% (~10%) 
for magnetic skyrmions in Pd/Fe/Ir (Pd/Pd/Fe/Ir) magnetic thin-hlms, which potentially 
could be detected in a CPP-device geometry. Establishing the physics of this generalized 
tunneling spin-mixing magnetoresistance (TXMR) could possibly inspire the design of future 
nano-magnetic devices based on such a mechanism. 


RESULTS 

System and procedure 

We study two magnetic thin-film heterostructures similar to Fig. lb purely from ab initio: 
fcc-overlayers of Pd/Fe and Pd/Pd/Fe on single-crystal bulk fcc-Ir(lll). These systems are 
attractive for a number of reasons. First, they generate large Dzyaloshinskii-Moriya interac¬ 
tions (DMIs) 36-38 , whose competition with the isotropic exchange interaction J determines 
the size and chirality of the skyrmions 4 . DMIs are large here because of the strength and 
nature of the inversion-symmetry-breaking in the heterostructures. At the Fe/Ir(lll) inter¬ 
face, a large spin-orbit interaction in the underlying heavy-metal substrate, here Ir(lll), is 
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relatively uncompensated by the overlayer Pd/Fe or Pd/Pd/Fe interface, leading to a large 
DMI vector preferentially in the plane of Fe, denoted by D. The ratio of |D|/J, along with 
an external magnetic field, can stabilize isolated skyrmions with diameter D$ k ~l-5 nm in 
size, and has been shown experimentally 20,23 . Second, by choosing a doublc-Pd overlayer 
(Pd/Pd/Fe/Ir) versus a single-Pd overlayer (Pd/Fe/Ir), one can alter the exchange interac¬ 
tions in Fe due to the modified nature of the interface hybridization and electronic charge 
transfer (see Supplementary Information SI). We investigate this effect to illuminate con¬ 
ceptual studies where other overlayer combinations and materials are used to engineer the 
size, shape, and stability of the isolated skyrmions 24,39,40 . 

We perform self-consistent density functional theory (DFT) calculations based on a full- 
potential Green function formalism including SOI 41,42 , which allows a perfect embedding of 
rea/-space defects, such as isolated skyrmions, into the ferromagnetic background system. 
Additional specifics of our computational scheme are given in the Methods section. 

Non-collinear inhomogeneity in nano-skyrmions 

Before coming to the essential physics of the tunneling spin-mixing magnetoresistance 
(TXMR) effect, we first self-consistently relax different-sized nano-skyrmions in otherwise 
ferromagnetic backgrounds (see Fig. 2), in both single- and doublc-Pd overlayer material 
stacks. We control the size of the skyrmionic defects by allowing different finite numbers 
of atoms to relax their magnetic moments in size and direction after the central atom has 
been spin-flipped as an initial condition. We investigate four different realistic skyrmion 
sizes: F>sk ~ 1-1, 1.7, 2.2, and 2.7 nm in diameter. The spin-textures exhibit a fixed and 
unique rotational sense as demanded by the DMI, which seeks energy gain by torquing 
the moments to rotate with respect to their neighbors. These structures are cycloidal and 
radial in nature as expected for magnetic thin-hlms. Thus our theoretical calculations are 
consistent in generating realistic nano-skyrmions which have been experimentally detected 
using magnetic spin-polarized currents 20,23 . 

We illustrate the spin-moment global rotation vs. the vertical (polar angle 6) of each 
atom and the pairwise difference between adjacent polar angles (d$). We will show that 
the spin-mixing perturbations to the local-density-of-states (LDOS) are a function of these 
angular parameters because the relative canting between different pairwise atomic sites 
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Figure 2 | Real-space relaxation of nano-skyrmions with increasing size, (a-c) Plots of 
axisymmetric cycloidal spin-whirls inside a magnetically active Fe-layer centered about 
increasingly larger skyrmionic defects in fcc-Pd/Fe overlayer on fcc-lr(lll) bulk substrate. 
Confining spins in the FM-background are shown transparent. We define 6 as the typical polar 
angle with the vertical and dd as the difference in polar angle between adjacent pairwise atoms, 
(d-f) Again but in fcc-Pd/Pd/Fe overlayer on fcc-lr(lll) bulk substrate. The smallest defects 
we study (D$k ~ 1-1 nm) are not shown here. The color bar in (a) represents for all 
illustrations. 


varies as a function of space, in addition to the absolute canting relative to the substrate. 
When analyzing each individual confined structure, we mention that d# itself is not constant 
between different pairs of atoms, such that there exists an inhomogeneity on the atomic- 
scale in the rotation of the magnetization direction with respect to the substrate plane. 
Furthermore, these inhomogeneities themselves are a function of diameter when comparing 
skyrmions of different sizes (Fig. 2a,b,c). 
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Electronic structure of isolated confined skyrmions 


We now move to establish the physics behind the TXMR effect within a scanning tunnel¬ 
ing microscopy/spectroscopy (STM or STS) experiment employing a non-spin-polarized tip, 
for which according to the Tersoff-Hamann model 43 , the differential conductance dl/dV is 
proportional to the LDOS of the sample, calculated at the tip position, R t i P , and the given 
bias energy E hias : 

^ocLDOS(R tip ;£ bia s,M) • (1) 

The LDOS depends on the configuration {s*} of atomic spins of the sample (i) relative to 
each other, e.g. in terms of d6 tJ = 6j — 6i for all atom-pairs (i,j), and (ii) relative to the 
lattice in terms of the absolute polar angle 0 t . The transport phenomenon related to the 
latter is known as the tunneling anisotropic magnetoresistance (TAMR) 33 , an effect related 
to spin-mixing due to the SOI. The former results from the spin-mixing hybridization of 
majority- and minority-states due to non-collinearity. Both can be subsumed as tunneling 
spin-mixing magnetoresistance (TXMR). Common to both is that the probability of tun¬ 
neling into majority- and minority-states depends on angles. The difference is that (i) is 
expected to be larger than (ii), since it is not a SOI effect. Also the appearance of both 
are different. For example, the TXMR due to (i) in a homogeneous magnetic spiral is the 
same across the spiral, because d6 VJ = d 9 for all atom-pairs (i, j), but different for spirals of 
different pitches under the transformation d 9 —> dd'. In contrast, the TAMR is modulated 
across the spiral 35 as 0* changes from atom-to-atom. Thus, the TXMR is used to measure 
conductance differences between two different magnetic states such as the difference between 
a skyrmion and the FM-state, but can also be used to resolve magnetization inhomogenieties 
inside complex spin-textures such as skyrmions or domain-walls. 

Since the spin-mixing perturbations due to non-collinearity and SOI are magnetic in 
nature, we show in Fig. 3a the spin-dependent LDOS in the magnetically-active Fe-layer as 
a function of the atomic position for the « 1.7 nm skrymion. For brevity we plot only 
the Pd/Fe/Ir case (see Supplementary Information S2 for the larger D $k ~ 2.2 nm case and 
S3 for the double-Pd cases). We note that the majority and minority spin-channels are given 
in the local spin-frame of reference for each Fe-atom. The color-coding of the plot, which 
corresponds to different atoms extending radially from the skyrmion’s center, is explained 
in Fig. 3c. The energy-zero is the Fermi energy, Ep = 0. 
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Figure 3 | LDOS in a realistic nano-skyrmion in Pd/Fe/lr. (a) LDOS in the 
magnetically-active Fe-layer resolved into minority and majority spin-channels. The 
resonance-peak near E ~ 0.5 eV in the FM-background (green) shifts in energy when 
approaching the center of the skyrmion (black), (b) The modification of the electronic 
structure in Fe contributes to a strong resonance in the LDOS in vacuum via hybridization 
through surface Pd-states. Arbitrary units are used so as to include in the same plot the nature 
of the Pd -d z 2 surface-state, whose resonance-peak near 0.5 eV survives in the vacuum. (c) 
Illustrative legend for (a, b) where the atoms are numbered and color-coded to identify the 
appropriate curves to the corresponding atoms in the skyrmion. The vacuum domains are 
represented by empty spheres. 

The resonant-states between 0.5 and 1.0 eV above E F are of Fe d-band minority-spin 
character and consist of d xy / x 2 _ y 2 , d xz / yz , and d z 2 states. These states hybridize with the 
sp states in the Pd-overlayer and give rise to Fe-Pd -spd z 2 hybrid states, named in short 
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Pd-d z 2 states localized in the Pd overlayer-film around E ~ 0.5 eV, as shown in Fig. 3b 
(black-dashed curve). It is clear that the surface-layer Pd-d z 2 state (shown only for the 
background-FM Pd-surface film), which has the proper orbital symmetry to decay slowly 
transverse to the substrate, controls the electronic structure in vacuum as a function of 
energy, characterized by a strong resonance in the vacuum-LDOS. An all-electrical STS 
measurement will be sensitive to this vacuum resonance-peak. 

The origin of this resonance and its behavior upon rotation can be understood by analyz¬ 
ing the energy window near E ~ 0.5 eV in Fig. 3a, where one can see the resonance-peaks 
shifting in energy in the Fe skyrmion-LDOS as a function of position. The green curve, 
which represents the ferromagnetic state of the background Fe-film, shows an electronic 
structure consistent with Fe-minority d-d hybridization when adjacent atoms couple ferro- 
magnetically (see Supplementary Section S2). Conversely, moving towards the center of 
the skyrmion, the quantization axes between two neighboring atoms becomes different, and 
majority-states of one atom can hybridize with minority-states of the second. This effect is 
especially pronounced at the central spin-flipped atom (black curve), where the resonance- 
peak has shifted lower in energy - as expected for antiferromagnetic (AFM) coupling (see 
Supplementary Section S2). We reproduce these effects within the context of a simple model, 
where we can qualitatively predict the change in LDOS as a function of the non-collinear 
magnetization rotation parameter dd as defined in Fig. 2a (see Supplementary Information 
S4 and S5). 

The energy-dependent disturbance to the LDOS resonance-peaks as a function of position 
moving radially along the skyrmion will manifest as a perturbation to the local electrical 
conductivity, and is the physical basis for the space-dependent TXMR effect. 

All-electrical skyrmion detection 

We now define the TXMR and make predictions for future experimental observation of 
the effect. The TXMR is the percent deviation of the local conductance from a reference 
conductance due to the spin-mixing from non-collinearity and SOI. As long as the magnetic 
state under consideration has a different non-collinearity than the reference state, there will 
be a TXMR. If one is interested in the spatial resolution of a complex spin texture (ignoring 
SOI), however, then an additional inhomogeneity within the non-collinearity is required, as 
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is the case for nano-skyrmions. 

The TXMR is by definition measured in vacuum. Here, we choose the reference to be 
somewhere far from the skyrmion in the FM-background. Then the normalized TXMR 
measured at site r is 


TXMR(r) 


LDOS^ - LDOS|g|(r) 

LDOS™ 


x 100% , 


( 2 ) 


where LDOSp^ is the LDOS in the vacuum just above the FM, and LDOS|| c j.(r) is the 
LDOS of the complex spin-texture in the vacuum just above site r. 

Integrating the TXMR over the entire device injection boundary, over all energies up to 
the bias energy eVbias, would give a measure of the total change in conductance, and would 
be the state-of-bit detection mechanism in a CPP-TXMR device like discussed in Fig. la. 
In an STS experiment, however, the effect could be amplified by selecting specific energy 
windows where the TXMR were largest as a function of position. 

In Fig. 4a,b we show the energy-resolved TXMR of the D$ k ~ 1-7 nm skyrmion’s central 
spin-flipped vacuum-site, with and without SOI, for the single- and doublc-Pd cases. We 
notice a sizeable TXMR effect for both systems. This holds true for all skyrmions that we 
studied, noting a small size-dependence of the effect which varies weakly as a function of 
D Sk (see Supplementary Information S6). 

TXMR signals in the different singlc-Pd and double-Pd material systems vary from 
skyrmion-to-skyrmion, however. Since the SOI arises from the Fe-Ir interface, the impact 
of SOI is much more pronounced in the single-Pd system (Fig. 4a), where the TXMR can 
peak out at an impressive ~40% when ignoring spin-orbit coupling, but decreases down to 
~ 20% when SOI is included. Interestingly, the spin-mixing due to SOI is to compete with 
the effects due to inhomogeneous non-collinearity, reducing the overall TXMR signal. In 
the double-Pd system, most of the TXMR signal is due to inhomogeneous non-collinearity, 
with a small contribution coming from the SOI (Fig. 4b). 

In Fig. 4c,d we plot the spacial variation of the TXMR signal, which differs significantly 
from atom-to-atom within the same skyrmion. We see that the TXMR effect is reduced 
when approaching the edge of the skyrmion (Fig. 4c,d blue curves), since the effective non- 
collinearity is reduced as the complex spin-texture fades into the ferromagnetic background. 

The vacuum-resonance we found in Fig. 3b appears now as a large TXMR signal at the 
same energy in the Pd/Fe/Ir(lll) system (Fig. 4a,c). Thus, an experimentalist probing the 
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Figure 4 | Tunneling spin-mixing magnetoresistance. (a,b) TXMR signals comparing the 
effects of SOI in a D Sk m 1.7 nm skyrmion, for single- and double-Pd systems. (c,d) Again for 
a line of atoms extending radially (see insets for labeling convention), (e) Expected STS-signal 
when approaching a skyrmionic defect in the single-Pd case. The electrical contrast has been 
projected onto the plane below the skyrmion. 
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surface in a STS experiment, having set the ripple-bias voltage near Vbias ~ 0.5V, would see 
an electrical contrast as visualized in Fig. 4e when approaching a skyrmion of similar size 
in the single-Pd heterostructure. 

Within a reasonable bias voltage range, the TXMR effect is smaller in the donble-Pd 
case (~10%) when compared to the single-Pd case (~20%). The additional Pd-overlayer 
changes the resonance nature around 0.5 eV compared to the single-Pd case, and states 
with a high tunneling cross-section into the vacuum are distributed over a wider energy (see 
Supplementary Information S3). As a consequence, the TXMR is reduced by nearly half 
when compared to the single-Pd case, and a bias energy near —0.5 eV is experimentally 
more favorable. 

DISCUSSION 

We have studied realistic- and experimentally-observable confined nano-skyrmions within 
metallic thin-hlms of Pd/Fe/Ir(lll) and Pd/Pd/Fe/Ir(lll) completely from first-principles. 
We established how the combined effects of local inhomogeneous magnetic non-collinearity 
and SOI in nano-skyrmions can alter the atomistic electronic structure in a magnetically 
active Fe-film, and, via hybridization with additional surface layers, the electrons which 
tunnel into the vacuum. 

The change in the LDOS can be understood in terms of the rotation parameters of 
the magnetic moment of the considered atom. The largest spin-mixing contribution comes 
from non-collinearity and depends on the relative canting between magnetic moments on 
neighboring sites, d 9. The dependence on the absolute polar angle of the magnetic moment 
with respect to the substrate, 9, comes in as a second-order term to the change in the LDOS, 
but can become important if the impact of the SOI is large. 

Finally, we have shown in detail how such a physical interplay could induce a sizeable elec¬ 
trical conduction anisotropy as a function of position and energy in realistic nano-skyrmions, 
up to ~20% in the single-Pd case. The manifestation of this TXMR effect could possibly 
be exploited in an all-electrical tunneling spectroscopy experiment. 

In addition, the changes in the magnetoresistance on the nanometer scale of skyrmions 
can possibly be engineered to design advanced magnetic memory devices. Typical memory 
circuits require at least one control device (either transistor or diode) in each memory cell. 
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Instead, technologies based on spin-mixing in single skyrmions could have potentially hun¬ 
dreds of bits stored in nanometer-sized magnetic racetracks needing only a single read-out 
element to detect the contents of each array (see Supplementary Information S7). 

METHODS 

The electronic structure was determined employing density functional theory (DFT) in 
the local spin density approximation 44 . Calculations were executed by means of the screened 
Korringa-Kohn-Rostoker full-potential scalar-relativistic Green function method 41,42 , where 
spin-orbit coupling was included self-consistently. A full-potential method is important to 
accurately describe the nature of the complex spin-texture and rapidly-decaying vacuum 
states of the tunneling electrons. 

For the calculations we chose an angular momentum cutoff of l max = 3 for the orbital 
expansions of the Green functions. The energy contour for numeric integration of the spin 
and charge density contained 40 grid points in the upper complex plane (including 7 Mat- 
subara poles) with a Brillouin zone mesh of 30x30 fc-points. The FM-slab LDOS and 
skyrmion impurity cluster LDOS were obtained by one-shot calculations using the FM-state 
or skyrmion-state as starting points, respectively. We found that increasing the 7-mesh to 
200 x200 was sufficiently adequate to numerically stabilize the relevant observables. 

The magnetic thin-klm slab configurations follow, where positive percentages refer to 
inward relaxation with respect to the Ir(lll) interlayer separation. We consider only fcc- 
stacking in all cases, which is in fact energetically favorable compared to hep crystal growth 39 . 

(a) Pd/Fe/Ir: 44 total layers (3 vacuum + 1 vacuum (-1%) + 1 Pd (8%) + 1 Fe (7%) + 
1 Ir (1%) + 33 Ir + 4 vacuum). 

(b) Pd/Pd/Fe/Ir: 44 total layers (3 vacuum + 1 Pd (-1%) + 1 Pd (8%) + 1 Fe (7%) + 
1 Ir (1%) + 33 Ir + 4 vacuum). 

We choose 34 Ir layers since it was the minimum thickness by which we completely 
decoupled any wave function penetration from top-to-bottom surface. We obtained the 
relaxation parameters as optimized and reported by Dupe et al 39 . 

To stabilize skyrmions after determining the 2D-FM slabs, the slab Green functions were 
harvested and a single spin-flipped Fe-atom was embedded in the FM-background. We 
then allowed three-layer cylindrical ring-like stacks of atoms within the skyrmion impurity 
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cluster to update their potentials and magnetic moments (Fe-layer + 1 Pd-layer above and 
1 Ir-layer below). The effect of the FM-background was included self-consistently by the 
slab Green function (Go), which connects the skyrmion impurity cluster (G imp ) to the host 
via the Dyson-like equation: Gi mp = Go + Go AV Gi mp , where AV represents the modified 
atomic potential as compared to the unperturbed slab Green function potential, V. In such a 
manner a real-space defect can be perfectly embedded in an otherwise periodic crystal. After 
converging the different sized skyrmionic profiles, observables were calculated as mentioned 
previously. 
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SUPPLEMENTARY INFORMATION 


SI — Charge density and angular momentum of FM thin-film heterostructures 


Below we show tables containing information regarding the charge density, charge trans¬ 
fer, and spin and orbital magnetic moments of the near-surface atomic layers in the fer¬ 
romagnetic (FM) state for each thin-film metallic system we studied, Pd/Fe/Ir(lll) and 
Pd/Pd/Fe/Ir(lll). The atomic number is Z, the total electronic charge is n, the charge 
transfer is An = (n — Z), the spin-moment is Mg, and the orbital-moment is Ml- Z, n, and 
An are given in units of electrons while Mg and Ml are given in units of /i-q. We give the 
layers from top-down in the singlc-Pd system. 


Pd/Fe/Ir(lll) 



Z 

n 

An 

Mg 

M l 

Pd 

46 

45.72 

-0.28 

0.31 

0.02 

Fe 

26 

26.18 

0.18 

2.70 

0.10 

Ir 

77 

76.84 

-0.16 

0.02 

-0.01 


We see that the spin magnetic moment of Fe is rather large (2.70 /ib) and the orbital moment 
is not negligible (0.10 hb)- Fe induces a sizable spin-moment in Pd of 0.31 hb which stiffens 
the isotropic exchange constant for first-nearest neighbors J compared to the system sl,S2 
having no Pd-overlayer (j Pd / Fc / Ir — 14.7 meV vs. J Fe / Ir = 5.7 meV). Pd is known to be 
a good Stoner system, i.e. it can easily develop a spin magnetic moment. Next we give a 
similar table for the double-Pd case (Pd/Pd/Fe/Ir): 


Pd/Pd/Fe/Ir(lll) 



Z 

n 

An 

Mg 

M l 

Pd 

46 

45.77 

-0.23 

0.08 

0.01 

Pd 

46 

45.94 

-0.06 

0.28 

0.02 

Fe 

26 

26.20 

0.20 

2.63 

0.09 

Ir 

77 

76.83 

-0.17 

0.02 

- 0.01 


Interestingly in this case, the spin magnetic moments of Fe and the nearest-neighbor Pd 
did not change greatly compared to the values obtained for Pd/Fe/Ir. However, the charge 
transfer of the nearest-neighbor Pd-layer is now much smaller. Since the surface Pd-layer 
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interacts with Fe indirectly through the inner Pd-layer, the induced spin magnetic moment 
is rather small (0.08 yU B ). 

S2 — LDOS in single-Pd 2.2 nm skyrmions 

Let us now turn to the case of a single-Pd system but for a larger skyrmion, Dm 2.2 nm, 
compared to the one shown in the main manuscript (i^sk ~ 1-7 nm, Fig. 3). In Fig. SI, we 
plot the LDOS in a realistic nano-skyrmion in Pd/Fe/Ir(lll) for the magnetically-active Fe- 
layer (a) and its corresponding vacuum sites (b), but with a slightly larger energy window 
[—1.5,+1.5] eV as compared to the main text, where a window of [—1.0,+1.0] eV was 
analyzed. 

Similar to the case of the smaller skyrmion in the single-Pd system as shown in the 
main manuscript (Fig. 3a, FM-Fe, green curve), a large resonance peak is observed in the 
background-FM Fe-hlm at 0.5 eV (Fig. Sla, FM-Fe, green curve). But by increasing the 
energy window beyond 1.0 eV, we clearly see a second resonance appearing at 1.25 eV in 
the nearly-FM atom Fe-3 at the edge of the skyrmion (Fig. Sla, Fe-3, navy curve). The 
appearance of these two resonance peaks in tandem suggests a strong interaction due to ferro¬ 
magnetic coupling. In other words, the resonant-peaks at 0.5 eV and 1.25 eV are actually the 
result of FM-interactions which have split a single resonance in two. An analogy can be made 
between these FM-splitted resonance-peaks and localized degenerate atomic orbitals at dif¬ 
ferent sites which upon interaction split in energy between bonding- and antibonding-states. 
The essential physics of this process can be reproduced within an Alcxander-Anderson model 
for d-d hybridization between FM-Fe d-states S3,S4 (see Supplementary Information S5). 

The splitting in energy between the two resonance-peaks reduces when the rotation angle 
of the magnetic moment increases (Fig. Sla, Fe-2, blue curve, and Fe-1, red curve). Fi¬ 
nally, in the central spin-flipped atom of the skyrmion (Fig. Sla, Fe-0, black curve), the 
splitting almost disappears as the two resonance-peaks quasi-merge. This is evidence of 
antiferromagnetic coupling between the central atom and its nearest-neighbors S3,S4 . We also 
reproduce this coupling within our Alexander-Anderson model in the antiferromagnetic case 
(see Supplementary Information S5). 

An important difference with respect to the center of the skyrmion (Fig. Sla, Fe-0, black 
curve) is observed when comparing to the smaller Dgk ~ 1-7 nm case (Fig. 3a, Fe-0, black 
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Figure S 1 | LDOS in a 2.2 nm skyrmion for Pd/Fe/lr(lll). a, With diameter Zl Sk « 
2.2 nm, shown is the LDOS in the Fe-layer resolved into minority (solid lines) and majority 
(dashed) spin-channels, b, LDOS in vacuum using arbitrary units so as to include the nature of 
the surface-layer Pd-d 2 2 state (dashed-black curve), c, The illustrative legend for (a, b) where 
the atoms are numbered and color-coded to identify the appropriate curves to the 
corresponding atoms in the skyrmion. The vacuum domains are represented by empty spheres. 

curve). In the smaller skyrmion, the quasi-merged resonance still shows two separate peaks, 
while in the larger skyrmion there is only a single, broad resonance peak. Thus we expect a 
difference in the TXMR signal in the larger skrymion compared to the D^k ~ 1.7 nm case. 
In general, there will be a weak size dependence of the spin-mixing signal, which should 
slowly decay as the diameter of the skyrmion is increased. In addition, from Fig. Sib we 
can see that the vacuum resonance near 0.5 eV in the smaller defect (Fig. 3b, solid curves) 
survives in the larger skyrmion (Fig. Sib, solid curves), and the TXMR remains detectible, 
as seen in Fig. S6a. 

To complete our analysis of the LDOS shown in Fig. SI, we mention that with regards 
to magnifying the strength of the TXMR effect, we need not restrict ourselves only to the 



LDOS in Fe 



LDOS in vacuum 
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energy window near 0.5 eV, but also look for other energy windows where a large change 
in the LDOS is observed as a function of position. We see this in the majority-states near 
—1.0 eV (Fig. Sla, dashed-curves), where the spin-mixing effect is perhaps even stronger 
than the 0.5 eV window composed of minority-states. Indeed, the TXMR signal is shown 
to be large for negative bias energies near —1.0 eV in this system (see Figs. 4a,c and S6a). 

S3 — LDOS of skyrmions in the double-Pd systems 

In Fig. S2, we plot the LDOS in a confined nano-skyrmion in Pd/Pd/Fe/Ir(lll) with 
diameter D$k ~ 1.7 nm for the magnetically-active Fe-layer (a) and its corresponding vacuum 
sites (b). 

Similar to the case of the smaller skyrmion in the singlc-Pd system as shown in the main 
manuscript (Fig. 3a, FM-Fe, green curve), a strong resonance-peak appears near 0.5 eV 
composed of minority-spin channel d-states (Fig. S2a, FM-Fe, green curve). This resonance- 
peak shifts in energy upon increased rotation of the magnetic moments moving towards the 
center of the skyrmion (Fig. S2a, Fe-2, blue curve, Fe-1, red curve, and finally Fe-0, black 
curve). This observed change in the electronic structure as a function of position due to the 
spin-mixing of majority- and minority-states suggests there will be a TXMR signal in this 
system as shown previously for the singlc-Pd system. 

Contrary to the singlc-Pd case, however, the resonance in the hybrid Fe-Pd-spd 2 2 state 
(Pd-d z 2 for short) observed in the Pd-overlayer around 0.5 eV vanishes, and instead appears 
near —0.5 eV (Fig. S2b, dashed-black curve). The steepness of this resonance leads to a 
disturbance in the vacuum-LDOS as a function of position in slightly lower energies, near 
—0.8 eV (Fig. S2b, solid curves). The character of the hybrid Pd-d 2 2 surface-state leads 
to a flat region in the vacuum-LDOS at positive energies. We then expect here a smaller 
TXMR signal in comparison to the single Pd-system for positive bias voltages. Therefore, 
in the double-Pd system, we suggest probing negative bias energies near —0.8 eV, where the 
spin-mixing effect is more significant, as exemplified in Fig. 4b,d. 

In Fig. S3, we plot the LDOS for a larger skyrmion, cs 2.2 nm, in Pd/Pd/Fe/Ir(lll). 
We use a similar labeling convention when decomposing the LDOS. Interestingly, the d- 
resonances are broader (Fig. S3a, solid curves) than for the smaller skyrmion (Fig. S2a, solid 
curves). In the center of the skyrmion, only one broad d-resonance is observed (Fig. S3a, 
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Figure S 2 | LDOS in a realistic 1.7 nm skyrmion in Pd/Pd/Fe/lr(lll). a, With 
diameter D Sk « 1.7 nm, shown is the LDOS in the Fe-layer resolved into minority (solid lines) 
and majority (dashed) spin-channels, b, LDOS in vacuum using arbitrary units so as to include 
the nature of the surface-layer Pd -d z 2 state (dashed-black curve), c, The illustrative legend for 
(a, b) where the atoms are numbered and color-coded to identify the appropriate curves to the 
corresponding atoms in the skyrmion. The vacuum domains are represented by empty spheres. 


Fe-0, black curve). To reiterate, this indicates that the TXMR signal will be different 
upon increasing the diameter of any skyrmion since its electronic structure is modified upon 
increasing D$k- From the shape of the vacuum-LDOS (Fig. S3b, solid curves), the TXMR 
is nicely detectable at an energy range around —0.8 eV, as shown in Fig. S6b, and in the 
same energy window as the smaller D$k ~ 1.7 nm Pd/Pd/Fe/Ir(lll) case. 

From Supplementary Sections S2 and S3 we can conclude that for each skyrmion in a 
Pd/Fe/Ir(lll) or Pd/Pd/Fe/Ir(lll) system, we can select an energy region in which a large 
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Figure S 3 | LDOS in a 2.2 nm skyrmion in Pd/Pd/Fe/lr(lll). a, With diameter 
.Dsk ~ 2.2 nm, shown is the LDOS in the Fe-layer resolved into minority (solid lines) and 
majority (dashed) spin-channels, b, LDOS in vacuum using arbitrary units so as to include the 
nature of the surface-layer Pd-d z 2 state, c, The illustrative legend for (a, b) where the atoms 
are numbered and color-coded to identify the appropriate curves to the corresponding atoms in 
the skyrmion. The vacuum domains are represented by empty spheres. 
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tunneling spin-mixing magnetoresistance (TXMR) will be found. Due to the relevant energy 
of detection, one can distinguish skyrmions of different sizes, and discriminate between 
skyrmions in Pd-films of different thicknesses. Finally, and perhaps most importantly, the 
spatial-variation of the spin-mixing signal means that the internal atomic structure of each 
individual skyrmion can be resolved. This means the TXMR can be used to visualize the 
size, shape, and structure of individaul defects. In devices based on spin-mixing, this also 
means that skyrmions can be used as nano-scopic information carriers, where the TXMR 
would be used to read the magnetic state-of-bit. 
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S4 — Angular dependence of the local-density-of-states in nano-skyrmions 

We now analytically derive the change in the LDOS at site r, denoted ALDOS(r), inside 
an axially symmetric skyrmion measured from the origin at r — 0 as a function of energy 
and magnetization rotation direction defined by the unit vector of the magnetic moment 
s(r) = ( sind(r) cos0(r), sin6 , (r) sin0(r), cos 0(r)). Examples of rotation parameters 9 and 
d 9 for given skyrmion magnetic moments can be found in Fig. S5d. When compared to the 
ferromagnetic state (with all moments pointing out-of-plane), there are two contributions to 
ALDOS(r), one due to spin-orbit interaction (SOI) and one due to non-collinearity (NC): 
ALDOS(r) = ALDOS so/ (r) + ALDOS Arc (r). The contribution from SOI, the so-called 
anisotropic magnetoresistance, is well known 85- . 

ALDOS so/ (r, E, 9) oc A(r, E) ■ [ 1 - s 2 z (r)] , (3) 

where A(r,E) is a coefficient depending on the site r and energy E, and s z (r) is the z- 
component of the spin-moment at site r. Thus upon including SOI, we expect, for example, 
a sin 2 d(r) dependence, which contributes to ALDOS(r) in second-order. 

The contribution from NC, intuitively, comes from the change in the electronic structure 
upon rotation of the magnetic moments at consecutive sites (i, j). For homogeneous mag¬ 
netic spirals, a constant deviation in the LDOS from the ferromagnetic (FM) state will be 
observed for each atom in the spiral (ignoring SOI). In such a spiral, the smooth rotation 
of moments 9{i) —» 9(j) = 9(i) + cl 9 for each atom pair is a symmetry operation commut¬ 
ing with a Hamiltonian having translational invariance, making each atom equivalent, and 
the electronic structure the same for each atom in the spiral. However, upon transforming 
the spiral’s pitch d 9 —» d 9', one would fold a different constant deviation of the electronic 
structure from the FM-state, such that spirals of different pitch can be identified by their 
different magnitudes in ALDOS^fo 

In skyrmions, however, the rotation of the magnetic moments is not homogeneous, i.e. 
d 9 is not constant for all atom pairs. Thus there will be a site-dependent deviation in the 
LDOS among the atoms inside a skyrmion with respect to each other. In what follows, we 
demonstrate that these deviations are a complex function of the rotation angles, which will 
depend on the details of the electronic structure, the energy probed, and even the size of 
the skyrmions. 
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To derive ALDOS^ (r), we utilize multiple scattering theory: G is the Green function 
describing the whole system upon rotation of the magnetic moments, and g is the Green 
function describing the initial FM-state. The Green function will be used to evaluate the 
change in the LDOS induced by the rotation of the magnetic moments: 

ALDOS iVC '(r, E) — —Tr L s [AG rr (E)] | , (4) 

as given in a matrix notation where a trace over orbital and spin angular momenta has to be 
performed. G can be evaluated via the Dyson equation connecting the non-collinear state 
to the ferromagnetic one: 

G — g + gAVG — g + gAVg + gAVgAVg + .... (5) 

where AV describes the change of the potential upon rotation of the magnetic moments. It 
can be expressed as: 

AD(r) = V dm (r) (a • s(r) - a z ) , (6) 

where <j is the vector of Pauli matrices, and Vdifj is the difference of the two spin components 
of the ferromagnetic potential (V ( ] — Vq)/2. 

We execute a similar expansion for the ferromagnetic initial Green function matrix g: 

g = <?suml 2 + gdiff&z , (7) 


where I 2 is the 2x2 identity matrix. 

Let us evaluate the first-order and second-order terms contributing to the Dyson Equ. 5: 

A G rr = g ri AV l g ir + g ri AV i g ij AV j g jr + .... (8) 

i ij 

where i and j are sites surrounding site r, or can be the site r itself. Since the trace over 
spin has to be performed, we will focus only on the terms that in the end will contribute to 
Equ. 4. We use a pair of useful properties of the Pauli matrices: 

Tr s [v x \ = Tr s [a y ] = Tr s [a z \ = 0 , (9) 

and 

(cr ■ s)(cr • s') = s ■ s' + i cr ■ (s x s') , (10) 
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( 11 ) 


where i is the imaginary unit. After simplifying, we find the following result: 

~A{Tr L s [g ri AV i g ir ] j = B rir (E) [1 - s*(i)] 
where the coefficient B is given by 

B rir = A{tYl [siVisS, + } . (12) 

In other words, the first sum in Equ. 8 leads to a behavior like (1 — cos#;). 

The second-order term is given by 

— ^{ T r L s \g ri AV i g ij AV j g jr ] } = C ri * [s(i) • § (j) - (s z (i) + s z (j)) + l] , (13) 

where the coefficient C is related to the Green functions and Vdiff by: 

C Hir = -l3(Tr L [j"„U' la jLVi9£ m + gSmVL&VLsi™ 

71 f 

+ alLvL&v^gZn + s2 ff vj f9 « m v3 e9 i] | . (14) 

Thus we obtain a dependence on the dot product of the unit vectors of the magnetic 
moments (1 — cosd#cosd0) and a contribution depending only on the ^-components of the 
unit vectors of the magnetic moments. 

We have thus demonstrated that due to NC, the dependence of the change in the LDOS 
with respect to the ferromagnetic state upon rotation of the magnetic moments is not trivial, 
and will have terms depending on the dot product between magnetic moments, contrary to 
the contribution coming from SOI. The non-collinear contribution is then 

ALDOS JVC '(r, E, {s}) = 

ZB rir (E)(l-s z (i))+ EC ri ^(E)[s{i)-s(j)-{s z {i) + s z (j)) + l] , (15) 

i %j 

where {s} is the spin configuration. Of course, depending on the details of the electronic 
structure and strength of perturbation related to the non-collinearity, higher-order terms 
can be important and have to be included in Equ. 15. 

Combining ALDOS^*" and ALDOS 507 , in the next section we will fit the change in 
the LDOS in terms of trigonometrical functions that depend on the rotation angles of the 
magnetic moments. We will apply these fits to our ab initio results as well as to an extended 
Alexander-Anderson model used to interpret the variation of the LDOS resonance-splitting 
upon rotation of the magnetic moments on neighboring sites. 
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S5 — Two-atom extended Alexander-Anderson model using Green functions 


We wish to estimate the change in the LDOS and qualitatively understand the shifting 
in energy of resonant d-states in Fe as a function of the rotation angle between adjacent 
moments. To this end, we consider for simplification two magnetic atoms (i, j ) = (1, 2) 
each having one localized orbital d z 2 whose single-particle eigenenergy is centered about 
E = e. The initial Hamiltonian describing this model is diagonal in spin-space. We could 
also consider an orbital of the type d xz in order to address the coupling induced by SOI 
between the d z 2 and d xz , as done by Caffrey et al ss . However, since the impact of SOI on 
the LDOS has already been discussed by others, we focus here on the impact of NO on the 
LDOS. We study the ALDOS as we vary dd = 9 \ — 62 between the two atoms. We restrict 
the hopping from atom-to-atom to non-spin-flip processes, characterized by the interaction 
parameter 14 op . 

In terms of Green functions, the following equation gives the LDOS for site 1: 

LDOS(l; E, {s}) = -^3 |tV s [G n (£)] j = jwg [E-H + iT}~ 1 j , (16) 

where T takes care of the broadening of the states. Instead of solving exactly the previous 
equation, one could also use perturbation theory, as described in the previous section S4, 
simplifying Equ. 15 to: 

ALDOS 7VC (l; E, {s}) = D(E) • (1 - cos d 9 ) , (17) 

where D = B 121 + C 1221 . 

The energy of the resonant d-states, their width, and splittings come from our first- 
principles calculations, e.g. Fig. Sla (atom Fe-3, navy curve). To obtain the proper splitting 
we choose a hopping parameter Vh op ~ 300-400 meV. We show the resulting LDOS in Fig. S4 
for five different rotation angles d 9 . There we reproduce the splitting of the resonance-peaks 
that we have seen in our first-principles calculations, where d-d hybridization is important, 
as seen in e.g. Fig. 3a of the main text, or Figs. Sla, S2a, and S3a of the Supplementary 
Information. 

Interestingly, for the extreme cases of a ferromagnetic state and an antiferromagnetic 
state, we recover the expectations of the Alexander-Anderson model S3,S4 . Indeed, when the 
two atoms are in the ferromagnetic state (dd = 0°), the d-d hybridization leads to a splitting 
of the original single orbital into bonding- and antibonding-states, seen as broad resonances 
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Figure S 4 | Alexander-Anderson prediction for d-d hybridization. Beginning from 
d#=0 o , the ferromagnetic state imposes a hybridization between the two-orbitals localized at 
the same energy. This produces a splitting into a bonding- and antibonding-state (green curve). 
Upon rotation, the atoms eventually become antiferromagnetically coupled (black curve). The 
change in the LDOS can be qualitatively estimated with a cos(dd) fitting parameter. 

near 0.5 eV and 1.25 eV, respectively (Fig. S4, green curve). The splitting is then given by 
2Vhop- If the magnetic state is antiferromagnetic (dd = 180°), there is repulsion between the 
minority and majority spin-states leading to a shift given by Vh 0 p/(Ef dz2 — Ej dz2 ). In our 
simple model, EE 2 — E. rJd2 is extremely large, thus the shift is not observed in Fig. S4. It 
is interesting to see how the splitting between the resonance-peaks decreases upon rotation 
of the magnetic moments until they merge to a single resonance in the antiferromagnetic 
case (Fig. S4, black curve). This is in accordance with the behavior of the LDOS calculated 
from first-principles for different sized skyrmions in the two systems Pd/Fe/Ir (Figs. 3a and 
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Sla) and Pd/Pd/Fe/Ir (Figs. S2a and S3a). 

Next, we wish to estimate the change in the LDOS as a function of rotation as previously 
discussed, which leads to the TXMR signal. In Fig. S5 we plot the change in LDOS compared 
against the background-FM for the model along with our ab initio results in Pd/Pd/Fe/Ir for 
skyrmions 1.7 nm and 2.2 nm in diameter, respectively. In addition we show fitted functions 
against 9 and c\6. Interestingly, a good fit to the change in LDOS shown in Fig. S5a,b is 
found by considering 

ALDOS(r; E, {s}) « A(E) ■ sin 2 9(r) + D(E) ■ [l-cos0(r)] (18) 

instead of strictly employing the terms given by Equ. 15. This result is similar to what 
we found in the simple two-orbital Alexander-Anderson model (Equ. 17). In such systems, 
the contribution from all spin-moment dot products (see Equ. 13) behaves on average like 
cos(0). This is naturally satisfied in the two-orbital model. 

In Fig. S5a the change in LDOS at the Fermi energy E = Ep = 0 is depicted. The 
blue dots represent the total ALDOS given from our ab initio calculations for each atom 
extending radially from the skyrmion core. The fitting-function f{9) is Equ. 18 and shown 
to be effective in fitting the first-principles data. We note that the change in the LDOS at 
the Fermi energy was not large. This is different than the situations in Fig. S5b,c where we 
probed the energy resonances near 0.5 eV, and a larger change in the LDOS was induced. 
In Fig. S5b, the fitting function f(6) is again Equ. 18 and shown to be effective in fitting 
the ab initio data- 

in Fig. S5c, we plot the ALDOS as given by the two-atom Alexander-Anderson model 
(black curve) with two fitting-functions. The first-order fit (red-dashed curve) based on 
Equ. 17 is shown to be slightly inaccurate when fitting the model ALDOS. Instead a 
higher-order term is needed to fit the data (blue-dashed curve). Thus we learned from 
the Alexander-Anderson model that depending on the probed energy, additional terms can 
be needed to improve the fit. A term proportional to sin 2 (d), similar to the one expected 
when SOI is included, improves considerably the fit to the non-collinear contribution to 
ALDOS. For example, near energy resonances, especially if they are sharp, the perturbative 
expansion up to second-order from Supplementary Section S4 begins to break down and 
higher-order terms are necessary to improve the fitting-functions. 
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(a) Pd/Pd/Fe/lr: D sk = 1.7 nm 
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(b) Pd/Pd/Fe/lr: D sk = 2.2 nm 


Figure S 5 | Energy-dependent change in LDOS in confined skyrmions. a, 

Atom-by-atom ab initio results (blue dots) of the total change in LDOS relative to the 
background-FM including SOI in 1.7 nm diameter skyrmion in Pd/Pd/Fe/lr fitted to f{9) 
(solid curve). The f(6) in (a,b) refers to Equ. 16, while the numbers labeling each dot 
respresent the corresponding atom in the inset extending radially from the skyrmion core, b, 
Again, first-principles results (blue dots) of the atomistic change in LDOS but for a 2.2 nm 
skyrmion in Pd/Pd/Fe/lr fitted against f(6). c, Fitting the Alexander-Anderson model for the 
non-collinear component of the ALDOS. d, Rotation parameters as found from the 
first-principles calculations, referring to (a) and (b). 


S6 — TXMR in larger 2.2 nm skyrmions 

Next we plot the TXMR signal spacially-resolved for the 2.2 nm skyrmions in both sys¬ 
tems. The spin-mixing signals retain the same general features and shapes as in the smaller 
defects (Fig. 4c,d). The important peak near 0.5 eV in the Pd/Fe/Ir case remains (Fig. S6a), 
along with another strong peak near —1.0 eV. In the case of Pd/Pd/Fe/lr (Fig. S6b), the 
TXMR has the strongest signal near —0.8 eV, as before. Therefore one conld infer from 
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Figure S 6 | Spacially-varying TXMR in 2.2 nm skyrmions. a, With diameter D S k ~ 2.2 
nm, we show the TXMR in vacuum for Pd/Fe/lr. The vacuum resonance near 0.5 eV 
(Fig. Sib) contributes a strong peak in the TXMR near the same energy, b, Again but in 
Pd/Pd/Fe/lr, where the peak-signal strength appears near —0.8 eV, as predicted by Fig. S3b. 

Fig. S6 that as the diameter of skyrmionic quasiparticles is increased, the spin-mixing effect 
not only survives, but also that specific locations of strong TXMR signals remain in similar 
energy windows as in the smaller structures. 

S7 — Skyrmion racetracks for dense magnetic memories 

Spin-transfer torque magnetic random access memory (STT-MRAM) circuits reliably 
read-out bit-states depending on a tunneling magnetoresistance anisotropy of ~30-50% in 
some structures 39 , with a hope to achieve a magnetoresistance ratio Ron/Roff ~ 200% 
by 2022 S1 °. A TXMR effect as large as r s-/ 20% as we have shown in this work should be 
enough to provide adequate read-margin for scaled technologies, and is larger than the < 2% 
change in resistance found in widespread commercially-used hard disk read heads based on 
anisotropic magnetoresistance alone sn . Smaller changes in magnetoresistance just means 
there should be a more sensitive read-out circuit. Typically this means a few extra control- 
and boost-transistors and does not substantially increase the footprint of the memory, i.e. 
incorporating more sensitive read hardware does not degrade packing density considerably. 

In potential skyrmion-based devices using CPP-TXMR, a -Rqn/-Roff ~ 120% could 
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Figure S 7 | Concept spin-mixing magneto-memory versus DRAM, a, 2x2 1-T 1-C 

DRAM unit. In a random access memory, any memory cell can be accessed in roughly the same 
time since the WLs and BLs are connected to each bit in parallel. The ideal minimum packing 
footprint goes as 4/ RRAM /bit in units of area, where /dram is the node generation, b, 
Spin-mixing racetracks are not random access but acquire a latency associated with moving the 
bits sequentially out of the way to reach the requested address. However, with realistic 
skyrmion velocities, the access time could be roughly the same as in DRAM, while dramatically 
increasing the packing density. 


feasably be well-worth the tradeoff when considering the possible performance gains with 
regards to: (1) potentially very low power dissipation due to small currents needed to ma¬ 
nipulate the magnetic textures; (2) fast speed operation due to reduced read/write latencies 
associated with nano-scopic size; and (3) large increases in packing density. Let’s consider 
points (2) and (3) in greater detail. 

One issue with racetrack memories is that they are not random access. In a random access 
memory (RAM), any read/write operation can access any bit with roughly the same access 
time since the word and bit access lines (WL and BL), which are connected to the set/reset 
elements and read elements, are also connected in parallel with the individual memory cells 
(see Fig. S7a). In a racetrack memory, the situation is different. In practice, if a read/write 
were requested for an address whose representative bit were at the end of the racetrack, 
there would be additional latency associated with moving subsequent magnetic domains out 
the way until the requested bit were under the read or write device (see Fig. S7b). But does 
this make racetrack memory intrinsically slower than RAM? 

This is not clear. It depends on the velocities of the skyrmions in the racetrack. Consider 
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200 nm length racetracks populated with 2 nm in diameter quasiparticles in an array of 100 
domains (or 100 bits). With velocity 312 t>sk = 100 m/s , the furthest bit would acquire an 
additional delay in time t delay 


fdelay 


Oil (til) (100 bits) 
100 (?) 


= 2 ns 


a small number compared to the total access time needed to complete the read or write 
operation, which in modern dynamic RAM (DRAM) is in the range 20-50 ns S13 . We do 
note, however, that larger in-plane currents will be required to accelerate the quasiparticles 
up to a velocity such as ugk = 100 m/s, meaning there will be a tradeoff between fdeiay and 
power consumption. 

Regardless, by incorporating skrymion racetracks based on spin-mixing, there seems to be 
at first glance negligible additional acquired access latency - in fact, we may learn to find in 
the end that racetracks can be potentially faster than traditional RAM in certain geometries 
and biasing conditions, due to the nano-scopic size of the skyrmion quasiparticles. 

With regards to circuit layouts, it is clear that moving to racetracks will provide large 
gains in packing density. As an example, let’s compare the workhouse 1-transistor 1- 
capacitor (1-T 1-C) DRAM unit to our racetrack-based spin-mixing magneto-memory. Con¬ 
sidering the lithographical node generation, or minimum feature size, /, DRAM memory 
minimum packing requirements for a single bit’s footprint in the ideal case goes as surface 
area SA DRA m ~ 4 f 2 per bit (Fig. S7a). The current technology node for DRAM in 2015 
is given as /dram ~ 22 nm by the International Technology Roadmap for Semiconductors 
(ITRS) S1 °. In a skyrmion racetrack, the potential effective per-bit SA could be reduced to 
possibly SArt ~ / RT P er bit, where /rt is the diameter of magnetic quasiparticles in the 
racetrack (Fig. S7b). In our study thus far we considered realistic skyrmions of order Dg k ~ 
2 nm. Comparing against current DRAM arrays, we define the potential gain in packing 


density V as 


*S74dram 4/5 ram (4) (22 nm) 


AA rt / 2 t (2 nm) 2 ' v ' 

Even including more intricate sense amplifiers for read/write operations in skyrmion race¬ 
tracks, a conservative estimate gives the increase in packing density around 10 x, with an 


upper limit around 500 x if incorporating ultimately-scaled sub-5 nm nanoskyrmions and 


common fabrication techniques. 
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Possibly the packing density could even be larger introducing vertical racetracks S14 , which 
are difficult to fabricate thus far. 
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